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Q ! Abstract 

Two-body charmless nonleptonic decays of the B s meson are studied within the framework of 
generalized factorization in which factorization is applied to the tree level matrix elements while 
the effective Wilson coefficients are p and renormalization scheme independent, and nonfactoriz- 
able effects are parametrized in terms of N^(LL) and N^ S (LR), the effective numbers of colors 
arising from (V — A)(V — A) and (V — A)(V + A) four-quark operators, respectively. Branching 
ratios of B s — > PP, PV, VV decays (P: pseudoscalar meson, V: vector meson) are calculated as a 
function of N^ S (LR) with two different considerations for N^ S (LL): (a) N^(LL) being hxed at 
the value of 2, and (b) Nf(LL) = Nf(LR). Tree and penguin transitions are classified into six 
different classes. We find that (i) the electroweak penguin contributions account for about 85% 



(for N^ S (LL) = 2) of the decay rates of B s — > rjir, rfir, rjp, rf p, 4>tt, 4>p, which receive contribu- 



tions only from tree and electroweak penguin diagrams; a measurement of them will provide a 
clean determination of the electroweak penguin coefficient ag, (ii) electroweak penguin corrections 
to B s — ► ur)(\ 4>ri,uj(j), <p,4>4> are i n general as significant as QCD penguin effects and even 
play a dominant role; their decay rates depend strongly on N° S (LR), (iii) the branching ratio of 
B s — ► f)ff i the analogue of B& — > t]'K, is of order 2 x 10 -5 , which is only slightly larger than that 
of rfrf, K* + p~ , K + K~ , K°K° decay modes, (iv) the contribution from the rf charm content is im- 
portant for B s — > rfrf, but less significant for B s — > r/r/' , and (v) the decay rates for the final states 
#+(*)#-(*) foH ow the pattern: T(B S K+K~) > T(B S -> K + K*~) > T(B S > 
T(B S -► K + *K-) and likewise for K^TC ( *\ as a consequence of various interference effects 
between the penguin amplitudes governed by the effective QCD penguin coefficients 04 and a^. 
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I. INTRODUCTION 



Recently there has been a remarkable progress in the study of exclusive charmless B de- 
cays, both experimentally and theoretically. On the experimental side, CLEO has discovered 



and a possible evidence for B — > <pK*. Moreover, CLEO has improved upper limits for many 
other channels. Therefore, it is a field whose time has finally arrived. On the theoretical 
aspect, many important issues have been studied in past years, such as the effective Wilson 
coefficients that are renormalization scale and scheme independent, nonfactorizable effects 
in hadronic matrix elements, the QCD anomaly effect in the matrix element of pseudoscalar 
densities, running light quark masses at the scale m&, and the q 2 dependence of form factors. 

In the present paper, we plan to extend previous studies of charmless hadronic decays of 
B~ , Bd mesons to the B s mesons. In principle, the physics for the B s two-body hadronic 
decays is very similar to that for the Bd meson except that the spectator d quark is replaced by 
the s quark. Experimentally, it is known that B ± — > r]'K ± and Bd — > v(K have abnormally 
large branching ratios, several times larger than previous predictions. It would be very 
interesting to see if the analogue of Bd — > tj'K, namely B s — > rji]' or B s — > rfrf still has 
the largest branching ratio in two-body B s charmless decays. Another point of interest is 
concerned with the electroweak penguin corrections. It is naively believed that in charmless B 
decays, the contributions from the electroweak penguin diagrams are negligible compared to 
the QCD penguin corrections because of smallness of electroweak penguin Wilson coefficients. 
As pointed out in [fj], some B s decay modes receive contributions only from the tree and 
electroweak penguin diagrams and moreover they are dominated by the latter. Therefore, 
electroweak penguins do play a dominant role in some of B s decays. There also exist several 
penguin-dominated B s decay modes in which electroweak penguin corrections to the decay 
rate are comparable to that of QCD penguin contributions. In this paper, we will study this 
in details. 

Experimentally, only upper limits on the branching ratios have been established for a few 
B s rare decay modes (see || or Table 7 of [||]) and most of them are far beyond the theoretical 
expectations. Nevertheless, it is conceivable that many of the B s charmless decays can be 
seen at the future hadron colliders with large b production. Theoretically, early systematical 
studies can be found in More recently, one of us (B.T.) || has analyzed the exclusive 

charmless B s decays involving the 7] or rf within the framework of generalized factorization. 

This paper is organized as follows. A calculational framework is set up in Sec. II in which 
we discuss the scale and scheme independent Wilson coefficient functions, parametrization 
of nonfactorizable effects, classification of factorizable amplitudes,..., etc. The numerical 
results and discussions are presented in Sec. III. Conclusions are summarized in Sec. IV. The 
factorizable amplitudes for all the charmless two-body B s decays are given in Appendices. 





II. CALCULATIONAL FRAMEWORK 
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A. Effective Hamiltonian 



The relevant effective AB = 1 weak Hamiltonian for hadronic charmless B decays is 

V '= \Vu b v: q ( Cl o? + c 2 o u 2 ) + v cb v: q ( Cl ot + c 2 o c 2 ) - v tb v* 1 

i=3 



H eS (AB = 1) = -L V uh V* uq {c x Ol + c 2 O u 2 ) + V A Vl{c x Ol + c 2 O c 2 ) - V tb V* cA + h.c, (2.1) 
V 2 i=3 

where q = d,s, and 



= (ub) v _ A (qu) v _ A , O u 2 = {qb) 

o 3 ( 5 ) = (qb) v _ A YX<H)v-aW+a)> °4( 6 ) = (q a b/3 



v _ A {uu) v _ A , 



v-a{v+a)i 



(2.2) 



°7(9) = 2 (l b )v-A £ e 9' (§V)v-m(v-a)> °8(10) 
9' 



^(q a b(3) v _ A J2 e i'Wa)v + A(v-A)^ 



with O3-O6 being the QCD penguin operators, O7-O10 the electroweak penguin operators, 
and (qiq2) v±A = 9i7/i(l ± 75)92- In order to ensure the renormalization-scale and -scheme 
independence for the physical amplitude, the matrix element of 4-quark operators has to 
be evaluated in the same renormalization scheme as that for Wilson coefficients q(/x) and 
renormalized at the same scale fi. Generically, the hadronic matrix element is related to the 
tree level one via 



<<?(//)) = <?(//) <0) tree , 



(2.3) 



with g(n) being the perturbative corrections to the four-quark operators renormalized at 
the scale \i. We employ the relation (2.3) to write (TC e s) = c eS (0) tree . Schematically, the 
effective Wilson coefficients are given by c eS = c([i)g(fi). Formally, one can show that c? 
are fi and renormalization scheme independent. It is at this stage that the factorization 
approximation is applied to the hadronic matrix elements of the operator O at the tree level. 
The physical amplitude obtained in this manner is guaranteed to be renormalization scheme 
and scale independent, f] 

Perturbative QCD and electroweak corrections to g(fi) from vertex diagrams and penguin 
diagrams have been calculated in |^-|TT|. The penguin-type corrections depend on k 2 , the 
gluon's momentum squared, so are the effective Wilson coefficient functions. To the next- 
to-leading order, we obtain [T2| 



„eff 



off 



„cff 



„cff 



.eff 



1.149, 

0.0211 + i0.0045, 
0.0134 + ^0.0045, 
-(0.0276 + i0.0369)a, 
-(1.318 + «0.0369)a, 



.off 



,,eff 



..eff 



off 



,,eff 
10 



-0.325, 

-0.0450 -iO.0136, 
-0.0560 - £0.0136, 
0.054 a, 
0.263 a, 



(2.4) 



*This formulation is different from the one advocated in @ in which the \i dependence of the 
Wilson coefficients c%(n) are assumed to be canceled out by that of the nonfactorization parameters 
Ss(n) and E\(n) so that the effective parameters af are \i independent. 
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at k 2 = ml/ 2. It is interesting to note that cf^ are very close to the leading order Wilson 
coefficients: c\° = 1.144 and c\° = -0.308 at fi= m b {m b ) |j| and that Re(cf_ 6 ) « |c^ 6 (/i). 



Therefore, the decay rates of charmless B decay modes dominated by QCD penguin diagrams 
will be too small by a factor of ~ (1-5) 2 = 2.3 if only leading-order penguin coefficients are 
employed for the calculation. 



B. Parametrization of nonfactorizable effects 

Because there is only one single form factor (or Lorentz scalar) involved in the class-I 
or class-II decay amplitude of B — > PP, PV decays (see Sec. II. C for the classification 
of factorizable amplitudes), the effects of nonfactorization can be lumped into the effective 
parameters ai and a 2 p 



af = of + cf (i- + xi) , af = cf + cf (±- + Xa ) , (2.5) 

where Xi are nonfactorizable terms and receive main contributions from color-octet current 
operators. Since | ff / c| ff | 3> 1, it is evident from Eq. (2.5) that even a small amount of 
nonfactorizable contributions will have a significant effect on the color-suppressed class-II 
amplitude. If Xi,2 are universal (i.e. process independent) in charm or bottom decays, then 
we have a generalized factorization scheme in which the decay amplitude is expressed in terms 
of factorizable contributions multiplied by the universal effective parameters af 2 . For B — ► 
VV decays, this new factorization implies that nonfactorizable terms contribute in equal 
weight to all partial wave amplitudes so that af 2 can be defined. It should be stressed that, 
contrary to the naive one, the improved factorization does incorporate nonfactorizable effects 
in a process independent form. For example, Xi — Xi — ~ \ i n the large- iV c approximation 
of factorization. Phenomenological analyses of the two-body decay data of D and B mesons 
indicate that while the generalized factorization hypothesis in general works reasonably well, 
the effective parameters af^ do show some variation from channel to channel, especially for 
the weak decays of charmed mesons |H-|TE||. An eminent feature emerged from the data 



analysis is that a| ff is negative in charm decay, whereas it becomes positive in the two-body 
decays of the B meson |l4] , p!7| ,|7| : 



• cS (D -> Kir) ~ -0.50, af (B -> Dtt) ~ 0.20 - 0.28. (2.6) 

It should be stressed that the magnitude of depends on the model results for form factors. 
It follows that 

X 2 {D -»• Ftt) ~ -0.36 , X2{B -»• -Dtt) ~ 0.12 -0.19. (2.7) 

The observation |x 2 (-B)| |x 2 (-D)| is consistent with the intuitive picture that soft gluon 
effects become stronger when the final-state particles move slower, allowing more time for 



significant final-state interactions after hadronization Phenomenologically, it is often to 
treat the number of colors N c as a free parameter to model the nonfactorizable contribution 
to hadronic matrix elements and its value can be extracted from the data of two-body 
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nonleptonic decays. Theoretically, this amounts to defining an effective number of colors 
Nf, called 1/f in 0, by 



1/Nf = {l/N c ) + X . (2.8) 

It is clear from (2.7) that 

Nf{D -> Kit) > 3, Nf(B -> Pvt) ~ 1.8 - 2.2 . (2.9) 

The effective Wilson coefficients appear in the factorizable decay amplitudes in the com- 
binations a 2 j = c|? + ^-c|f_ 1 and a 2 j_i = c|f_ 1 + ^-c|f (z = l,---,5). As discussed in 
the Introduction, nonfactorizable effects in the decay amplitudes of B — > PP, VP can be 
absorbed into the parameters a| ff . This amounts to replacing N c in a| ff by (Nf)i. Explicitly, 

a 2f = C 2f + / 7\r P fn C 2i-1' °2i-l = C 2i-1 + / Ar P ff\ C 2i ) (* = 1) ' ' ' > ^). (2.10) 

It is customary to assume in the literature that (Nf) 1 ~ (Nf) 2 • - • ~ (iV r ^ ff ) 10 so that the 
subscript z can be dropped; that is, the nonfactorizable term is usually assumed to behave in 
the same way in penguin and tree decay amplitudes. A closer investigation shows that this 
is not the case. We have argued in Jl2[ that nonfactorizable effects in the matrix elements 
of (V — A)(V + A) operators are a priori different from that of (V — A)(V — A) operators. 
One reason is that the Fierz transformation of the (V — A)(V + A) operators 05,6,7,8 is quite 
different from that of (V — A)(V — A) operators Oi, 2,3,4 and O 9 ,i . As a result, contrary to the 
common assumption, Nf(LR) induced by the (V — A)(V + A) operators are theoretically 
different from Nf(LL) generated by the (V — A)(V — A) operators [jl^] . From Eq. (2.10) it 
is expected that 

Nf(LL) ee (Nf) i » (Nf) 2 » (Nf ) 3 » (iVf ) 4 » (iVf ) g » (iVf ) m , 
Nf(LR) = (Nf) 5 « (iVf ) 6 » (iVf ) 7 « (iVf ) g , (2.11) 

and Nf(LR) ^ Nf(LL) in general. In principle, Nf can vary from channel to channel, as 
in the case of charm decay. However, in the energetic two-body B decays, Nf is expected 
to be process insensitive as supported by data J1J. 

The N c eff -depend ence of the effective parameters af's are shown in Table I for several 
representative values of Nf. From Table I we see that (i) the dominant coefficients are 
ai, a 2 for current-current amplitudes, 04 and for QCD penguin- induced amplitudes, and 
ag for electroweak penguin- induced amplitudes, and (ii) ai, a 4 , and ag are iV^-stable, while 
others depend strongly on Nf. Therefore, for charmless B decays whose decay amplitudes 
depend dominantly on iV^-stable coefficients, their decay rates can be reliably predicted 
within the factorization approach even in the absence of information on nonfactorizable 
effects. 

The CLEO data of B^ — ► ojtx^ available last year clearly indicate that Nf(LL) is favored 
to be small, Nf(LL) < 2.9 |12|. If the value of Nf(LL) is fixed to be 2, the branching 
ratio of B ± — > um^ for positive p (p being a Wolfenstein parameter; see Sec. II. D), which is 
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preferred by the current analysis [ID ], will be of order (0.9 — 1.0) x 10~ 5 , which is very close 
to the central value of the measured one. Unfortunately, the significance of B ± — > ujtt^ 1 is 
reduced in the recent CLEO analysis and only an upper limit is quoted [pD[ . Nevertheless, 
the central value of BiB^ — > 7T o>) remains about the same. Therefore, a measurement of 
its branching ratio is urgently needed. A very recent CLEO analysis of B° — > 7r + 7r~ plf 
presents an improved upper limit, B(B° — > 7r + 7r~) < 0.84 x 10~ 5 . If the form factor P^P^O) is 
known, this tree-dominated decay could offer a useful constraint on Nfr(LL) as its branching 
ratio increases slightly with Nf. For Ff*(0) = 0.30, we find Nf(LL) < 2.0. The fact that 
N^ S (LL) is favored to be at the value of 2 in hadronic charmless two-body decays of the 
B meson is consistent with the nonfactorizable term extracted from B — > (D,D*)tt, Dp 
decays, namely N^ S (B — > Dtt) « 2. Since the energy release in the energetic two-body 
decays B — > um, B — > Dti is of the same order of magnitude, it is thus expected that 
N^ s (LL)\b^ujtt ~ 2. In analogue to the class-Ill B — > Dtt decays, the interference effect of 
spectator amplitudes in charged B decays B~ — > 7r~7r°, p~n°, n~ p° is sensitive to N^ S (LL); 
measurements of them will be very useful to pin down the value of (LL). 

As for N^ G (LR), it is found in |E| that the constraints on N^ S (LR) derived from B ± — > 



(j)K and B — > <pK* are not consistent. Under the factorization hypothesis, the decays 
B — > <pK and B — > <pK* should have almost the same branching ratios, a prediction not 
borne out by current data. Therefore, it is crucial to measure the charged and neutral 
decay modes of B — > <p(K,K*) in order to see if the generalized factorization approach is 
applicable to B — ► cf)K* decay. Nevertheless, the analysis of B — > rj'K in [|TJ] indicates that 
N° G (LL) w 2 is favored and N* S (LR) is preferred to be larger. Since the energy release in 
the energetic two-body charmless B decays is not less than that in B — ► Dtt decays, it is 
thus expected that 

|x(2 - body rare B decay) | < \x(B -> Dir)\. (2.12) 

It follows from Eqs. (2.7) and (2.8) that Nf(LL) w Nf(B -> Dn) ~ 2 and Nf(LR) ~ 2-5, 
depending on the sign of \- Since N* (LR) > N° G (LL) implied by the data, therefore, we 
conjecture that 

Nf(LL) w 2, iVf < 5. (2.13) 



C. Factorizable amplitudes and their classification 

Applying the effective Hamiltonian (2.1), the factorizable decay amplitudes of B s — ► 
PP, VP, VV obtained within the generalized factorization approach are summarized in the 
Appendices A,B,C, where, for simplicity, we have neglected P^-annihilation, space-like pen- 
guins and final-state interactions. All the penguin contributions to the decay amplitudes can 
be derived from Table II by studying the underlying b quark weak transitions. To illustrate 
this, let x( BMl,M2 > denote the factorizable amplitude with the meson M 2 being factored out: 

X (BM U M 2 ) = ( M 2 |(g 2 g 3 )y-A|0)(M 1 |(g 1 6) v ,_jP). (2.14) 
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Table I. Numerical values for the effective coefficients af at Nf = 2, 3, 5, oo (in units of 10 4 for 
■ ■ ■ , aio)- For simplicity we will drop the superscript "eff" henceforth. 





iV c cfr = 2 


Nf = 3 


iV c cff = 5 


Nf = oo 


ai 


0.986 


1.04 


1.08 


1.15 


a 2 


0.25 


0.058 


-0.095 


-0.325 


a 3 


-13.9 - 22.6i 


61 


121 + 18.1?' 


211 + 45.3?' 


0>4 


-344 - 113i 


-380 - 121?' 


-408 - 127?' 


-450 - 136i 


a 5 


-146 - 22.6i 


-52.7 


22.0 + 18.1? 


134 + 45.3? 


a 6 


-493 - 113i 


-515 - 121? 


-533 - 127i 


-560 - 136? 


a 7 


0.04 — 2.73? 


(\ rj -i i) TO ' 

— 0.71 — 2.73? 


-i o A O TO ' 

— 1.24 — 1.161 


—2.04 — 1.161 


a 8 


2.98 - 1.37? 


3.32 - 0.9H 


3.59- 0.55?' 


4 


a 9 


-87.9 - 2.73?' 


-91.1 - 2.73?' 


-93.7- 2.73?' 


-97.6- 2.73i 


aio 


-29.3 - 1.37? 


-13.1-0.91? 


-0.04 - 0.55? 


19.48 



In general, when M 2 is a charged state, only a even penguin terms contribute. For example, 
from Table II we obtain 

A(B S -> ^ + 7r-) pcng oc [a 4 + an, + (a 6 + a 8 )R] X^ BsK+ ^~\ 
A(B S -> X* + 7r-) peng oc [a 4 + a 10 - (a 6 + a 8 )#] x< B '** + ' w ~>, 
A(5 S -> tfVW oc [a 4 + aio] X^ + '"~), (2.15) 

with R' ^ R ^ m%/ (mbrrid) ■ When M 2 is a neutral meson with I 3 = 0, namely, M 2 = 7r°, p°, a; 
and a Q dd penguin terms start to contribute. From Table II we see that the decay 
amplitudes of ~B S -> Mn°, ~B S -> Mp°, U s -> Mcj, U s -> M77O contain the following 
respective factorizable terms: 

\{-a 7 + a 9 )X^\ 
\{a 7 + a 9 )X^\ 
(2a 3 + 2a 5 + 1 -a 7 + ^a 9 )Xi B ^\ 
{2a,-2a,-\a 7 + \a 9 )X^ M ^' ) \ 



where the subscript u indicates the uu quark content of the neutral meson: 

X (B S M,„°) = (7r O| Wy _ A | 0>(Mi | (? - i6)v _J^ ) _ 



For example, the penguin amplitudes of B s — > 77a; and X n are given by 



V- 

0^0 



(2.16) 
(2.17) 



A(B S -> r/^Opcng oc 



2(a 3 + a 5 ) + i(a 7 + a 9 ) 



A(5 S - #V) peng oc |(-a 7 + a 9 )X^ K °^ + 



y(B s ri,w) 

1 . 1 , ■ 

a 4 - -aio + (a6 - -a 8 )-R 



OC 



3 1 1 



(2.18) 
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respectively. It is interesting to note that the decays B s — > (rf\ 0)(vr°, p°) do not receive any 
contributions from QCD penguin diagrams and they are dominated by electroweak penguins. 
We will come back to this interesting observation later. 

Table II. Penguin contributions to the factorizable B — > PP, VP, VV decay amplitudes multiplied 
by -{G F /V2)V tb V t *, where q = d,s. The notation B -> M u M 2 means that the meson M 2 can be 
factored out under the factorizable approximation. In addition to the a even terms, the decay also 
receives contributions from a Q( jd penguin effects when M 2 is a neutral meson with I3 = 0. Except 
for 77 or rf production, the coefficients R and R' are given by R = 2mp/[(mi + m2)(m;, — 777-3)] and 
R' = — 2m? P /\{m\ + m 2 )(mb + 777,3)], respectively. 



Decay 


b — > aim, 6 — ► gcc 


b — > gdd, 6 — > gss 


B -> P,P 


04 + aio + («6 + G8)-R 


04 — ^aio + («6 — ^o,s)R 


B -^V,P 


04 + aio + («6 + as)R' 


04 — ^aio + (a6 — 2 a s)R' 


B —> P,V 


a 4 + aio 


(14 — ^aio 


B -> V,V 


a 4 + aio 




B -» P,P° 


03 - a 5 - a 7 + a 9 


«3 - «5 + ^«7 - ^«9 


B ^V,P° 


03 - a 5 - a 7 + a 9 


03 - a 5 + |a 7 - ±a 9 


B -> P,V° 


03 + 05 + a 7 + ag 


^3 + «5 - |a 7 - ^«9 


B -> V, V° 


03 + «5 + 07 + «9 


«3 + «5 - \a 7 - ±a 9 



Just as the charm decays or B decays into the charmed meson, the tree-dominated 
amplitudes for hadronic charmless B decays are customarily classified into three classes [|THJ : 

• Class-I for the decay modes dominated by the external P^-emission characterized by 
the parameter a\. Examples are B s — > K + n~ , K* + n~, ■ ■ -. 

• Class-II for the decay modes dominated by the color-suppressed internal P^-emission 
characterized by the parameter a 2 . Examples are B s — > K°ir , K°p°, 

• Class-III decays involving both external and internal W emissions. Hence the class-III 
amplitude is of the form a\ + ra 2 . This class does not exist for the B s . 

Likewise, penguin-dominated charmless B s decays can be classified into three categories: f\ 

• Class- IV for those decays whose amplitudes are governed by the QCD penguin param- 
eters 04 and ci6 in the combination 04 + Ra^, where the coefficient R arises from the 
(S — P)(S + P) part of the operator Oq. In general, R = 2m P /[(mi + m 2 ){mb — m 3 )] 
for B — ► P a Pb with the meson being factored out under the factorizable ap- 
proximation, R = —2mp /[(mi + m 2 )(mj + 777,3)] for B — > V a Pb, and R = for 



^Our classification of factorizable penguin amplitudes is not the same as that in [|22|; we introduce 
three new classes in the same spirit as the classification of tree-dominated decays. 



B — > P a H and 5 — > V^Vf,. Note that is always accompanied by aio, and a6 by 
a§. In short, class-IV modes are governed by a eve n penguin terms. Examples are 
B s -> if+if-, 0r/O, • • 

• Class-V modes for those decays whose amplitudes are governed by the effective coef- 
ficients 0,3, 0,5, a 7 and a 9 (i.e. a Q( jd penguin terms) in the combinations a 3 ± a 5 and/or 
a-? ± a$ (see Table II). Examples are B s — > 7rr/'), wr/), 7r0, • • •. 

• Class- VI involving the interference of class-IV and class-V decays, e.g. B s — ► 

Sometimes the tree and penguin contributions are comparable. In this case, the inter- 
ference between penguin and spectator amplitudes is at work. There are three such decays: 
B s — > K°lu, K* Q rf\ K*°lu; they involve class-II and -VI amplitudes (see Tables IV and V). 

D. Input parameters 

In this subsection we specify the values for various parameters employed in the present 
paper. For current quark masses, we employ the running masses at the scale \i = m b . 

m u (m b ) = 3.2 MeV, m d (m b ) = 6.4 MeV, m s (m b ) = 105 MeV, 

m c (m b ) = 0.95 GeV, m b (m b ) = 4.34 GeV (2.19) 

As for the Wolfenstein parameters A, A, p and 77, which are utilized to parametrize the quark 
mixing matrix, we use A = 0.804, A = 0.22, p = 0.16 and 77 = 0.34. The values for p 
and 77 follow from a recent analysis of all available experimental constraints imposed on the 
Wolfenstein parameters |fl9[] : 

p= 0.156 ±0.090, 7/ = 0.328 ± 0.054, (2.20) 

where p = p(l — 4r) and fj = 7/(1 — 4r). For the values of decay constants, we use f n = 132 
MeV, f K = 160 MeV, f p = 210 MeV, f K * = 221 MeV, f u = 195 MeV and ^ = 237 MeV. 

To determine the decay constant f^ { , v defined by (0 1?7^75?| rf^) = if q ^)Pij,, it has been 
emphasized [0,0 that the decay constants do not simply follow the 77 — 7/ state mixing 
given by 

rf = 7/s sin 9 + 770 cos 9, 77 = rjs cos 9 — r] sin 9. (2-21) 

Introduce the decay constants /§ and / by 

(OK|77o> = i/oP M , (0\Al\ Vs ) = if sPll . (2.22) 

Because of SU(3) breaking, the matrix elements (0|A^ 8 ^ 1 778(0) ) do not vanish in general and 
they will induce a two-angle mixing among the decay constants, that is, f% and fL are 
related to / 8 and f by 



9 



Likewise, 



^ = VE sin 08 + {k cos 6>0 ' 



n = — cos 6* 8 - sin 00, 



^ = ~ 2 vl sin 08 + {h cos 6 °' 



^ = ~ 2 cos ° 8 ~ vl sin 00 ' 



(2.23) 



(2.24) 



Based on the ansatz that the decay constants in the quark flavor basis follow the pattern of 
particle state mixing, relations between 9s, #o and 9 are derived in ||25|| , where 9 is the T] — r]' 



mixing angle introduced in (2.21). It is found in |25| that phenomenologically 
9 8 = -21.2°, 9 = -9.2°, 9 = -15.4°, 



and 



is/A = 1.26, /„/% = 1.17. 



(2.25) 



(2.26) 



The decay constant fZ, defined by (0|c7 M 75c|?/) = if^q^, has been determined from 
theoretical calculations p7|-p9|1 and from the phenomenological analysis of the data of 
J /if) — > rjcj, J /if) — > 7/7 and of the 777 and 7/7 transition form factors [11,25,30 



it 



lies in the range -2.3 MeV < < -18.4 MeV. In this paper we use the values 
= -(6.3 ± 0.6) MeV, /° = -(2.4 ± 0.2) MeV, 



(2.27) 



as obtained in 25 



For form factors, the Bauer-Stech-Wirbel (BSW) model [26] gives |5j fj] 



F B * K (0) 


= 0.274, 


F Bs ^(0) 


= 0.335, 


F<f^(0) 


= 0.282, 




Ag'+{Q) 


= 0.272, 




= 0.273, 


A?'*(0) = 


= 0.273, 






= 0.236, 


A B ° K *(0) 


= 0.232, 




= 0.231, 






= 0.319, 


V B ' K '(Q) 


= 0.281. 






(2.28) 



It should be stressed that the 77 — 77' wave function normalization has not been included in the 
form factors Fq" 1 " and F /as ~; they are calculated in a relativistic quark model by putting 
the ss constitutent quark mass only. To compute the physical form factors, one has to take 
into account the wave function normalizations of the rj and rj': 



^0 



,V6 



cos 9 



sin 9 Fq 



B s r) S s 



iB a rf 



^__L sin ^ + _L cos ^ F( f^. (2.29) 



It is clear that the form factors F BsV and Fq" 11 have opposite signs. 



J The form factors adopted in are calculated using the light-front quark model and in general 
they are larger than the BSW model's results. 
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For the q 2 dependence of form factors in the region where q 2 is not too large, we shall 



use the pole dominance ansatz, namely, 



/(0) 



, -i ■> o\n i {2.60) 

(1 - q 2 /m 2 ) 

where m* is the pole mass given in |18|| . A direct calculation of B —>■ P and B — > V form 
factors at time-like momentum transfers is available in the relativistic light-front quark model 
with the results that the q 2 dependence of the form factors A , A 2 , V, Fi is a dipole 
behavior (i.e. n = 2), while Fq, A\ exhibit a monopole dependence (n = 1). 

Recently, the B s — > K* and B s — > <f> form factors have also been calculated in the 
light-cone sum rule approach [^4] with the parametrization 

/(0) 



1 — a(q 2 /m 2 Bs ) + b(q 2 /m 



(2.31) 



for the form-factor q 2 dependence. The results are 



A B ^(0) 


= 0.382, 


a 


= 1.77, 


b = 


0.856, 


A Bs (O) 


= 0.296, 


a 


= 0.87, 


b = 


-0.061 


Af s0 (O) 


= 0.255, 


a 


= 1.55, 


b = 


0.513, 


V B ^(0) 


= 0.433, 


a 


= 1.75, 


b = 


0.736, 


A** K *(0) 


= 0.254, 


a 


= 1.87, 


b = 


0.887, 


A^ K *(0) 


= 0.190, 


a 


= 1.02, 


b = 


-0.037 


A*° K *(0) 


= 0.164, 


a 


= 1.77, 


b = 


0.729, 


V B ° K \0) 


= 0.262, 


a 


= 1.89, 


b = 


0.846 . 



(2.32) 

It is obvious that the q 2 dependence for the form factors A , A 2 and V is dominated by 
the dipole terms, while Ai by the monopole term in the region where q 2 is not too large. 
In Tables IV and V we will present results using these two different parametrizations for 
B s — > V form factors. 

We will encounter matrix elements of pseudoscalar densities when evaluating the penguin 
amplitudes. Care must be taken to consider the pseudoscalar matrix element for r/'- 1 — ► 
vacuum transition: The anomaly effects must be included in order to ensure a correct chiral 



behavior for the pseudoscalar matrix element ||12|| . The results are |35| . p 

2 



2m« 



(t/'>|u7 5 u|0) = (v {,) \d-f 5 d\0) = r (0 (rj^\s l5 s\0), 



(2.33) 



with W2 



2f 2 -f 2 cos9 + ± S m9 



2/1 



2/ 2 



/2 cos6» - v^sin^' 
fi cos 6 — y/2 sin 9 



2/| -/o 2 cos^ + ^sin^ 



(2.34) 
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III. NUMERICAL RESULTS AND DISCUSSIONS 



With the factorizable decay amplitudes summarized in Appendices and the input pa- 
rameters shown in Sec. II, we are ready to compute the branching ratios for the two-body 
charmless nonleptonic decays of the B s meson. The decay rates for B s — > PP, VP are given 
by 

T(B, -> VP) = -^-\A(B, -> VP)l(e -p B )\\ (3.1) 



57rm 



The decay B s — > VV is more complicated as its amplitude involves three form factors. In 
general, the factorizable amplitude of B s — > V{V 2 is of the form: 

A (B S - V X V 2 ) = aXPMM + (3X^ V ^ 

= ( ai A?'* + fhAf-^el ■ e* 2 + (a 2 A B ^ + f3 2 A^)(e{ ■ pj(e* 2 ■ pj 
+ ie^sTsTp p Ba P a i(^ vBBVl + foV B ° v >), (3.2) 

where use of Eq. (CI) has been made. Then 

T(B S - V X V 2 ) = -^\a 1 (m Ba +m 1 )m 2 f V2 A B ^(m 2 2 )\ 2 (H + 2(H 1 + 2( 2 H 2 ), (3.3) 

where 

H = (a-bxf + 2(1 + cV), 
H\ = (a — bx)(a — b'x) + 2(1 + cc'yy'), 

H 2 = (a-b'x') 2 + 2(l + c ,2 y' 2 ), (3.4) 

with 

m B s -m 2 -ml 2m 2 Bs p 2 c 2m Ba p c 



2m x m 2 ™#2( m B s + m i) 2 ' { m B s + m-i) 2 ' 



Tfln 



C a x A BaVl (ml) ' * Af^(mi)' V Af'*(mg)' (3 ' 5) 

where p c is the cm. momentum, m x (m 2 ) is the mass of the vector meson V x [V 2 ), and 
6', c', x', y' can be obtained from b, c, x, y, respectively, with the replacement V x <-> V2. 

The calculated branching ratios for 5 S . — ► PP, VP, VV decays averaged over CP- 
conjugate modes are shown in Tables III-V, respectively, where the nonfactorizable effects 
are treated in two different cases: (i) N^ S (LL) 7^ N^ S (LR) with the former being fixed at the 
value of 2, and (ii) Nf(LL) = Nf(LR). For decay modes involving B s -> K* or B s -> 
transition, we apply two different models for form factors: the BSW model [see (2.28)] and 
the light-cone sum rule approach [see (2.32)]. To compute the branching ratio, we have used 
the B, lifetime [Hi 
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t(B s ) = (1.54 ± 0.07) x 10- u s. (3.6) 

From Tables III-V we see that the branching ratios for class-I and -IV modes are stable 
against the variation of as they depend on the coefficients a\, and a§ which are N* - 
stable. Class- V channels in general depend on the coefficients 03 + a 5 and a 7 + a 9 . However, 
the decays 

B s -> T]TT, 77V, T]p, Tf'p, 07T, 0p (3.7) 

do not receive any QCD penguin contributions 0. Therefore, these six decay modes are 
predominantly governed by the electroweak penguin coefficient ag, which is A^ ff -insensitive. 
A measurement of them can be utilized to fix the parameter ag. Note that their branching 
ratios are in general small, ranging from 4 x 10~ 8 to 0.4 x 10~ 6 , but they could be accessible 
at the future hadron colliders with large b production. 

In order to see the relative importance of electroweak penguin effects in penguin- 
dominated B s decays, we follow p2l to compute the ratio 



B(B S ^> h 1 h 2 )(vnth a 7 , ■■■ J a w = 0) 

Rw ~ b( Bs =; hm • (3 - 8) 

Obviously, if the tree, QCD penguin and electroweak penguin amplitudes are of the same sign, 
then (1 — R w ) measures the fraction of non-electroweak penguin contributions to B(B S — > 
hih 2 )- It is evident from Table VI that the decays listed in (3.7) all have the same N* s - 
dependence: For N° S (LL) = 2, the electroweak penguin contributions account for 85% of 
the branching ratios for B s —>■ i]7r, ■ ■ ■ , 0p, and the ratio Rw is very sensitive to Nfr when 
N^i^LL) = N° S (LR). We also see that electroweak penguin corrections to 

B s — > u>r], ur]', <pr], 0r/, u>(f), K<fi, K*<fi, 00, (3.9) 

depending very sensitively on N^ s , are in general as important as QCD penguin effects and 
even play a dominant role. For example, about 50% of B(B S — > K <p) comes from the 
electroweak penguin contributions at Nf(LL) = 2 and Nf(LR) = 5. 

Strictly speaking, because of variously possible interference of the electroweak penguin 
amplitude with the tree and QCD penguin contributions, Rw is not the most suitable quan- 
tity for measuring the relative importance of electroweak penguin effects. For example, it 
appears at the first sight that only 21% of B(B S — > wq') and B(B S — > ui(j>) arises from the 
electroweak penguins at N° a (LL) = N^ S (LR) = 3. However, the decay amplitudes are 
proportional to (see Appendix B) 

V ub V* s a 2 - V tb V*[2(a 3 + a 5 ) + ^(a 7 + a 9 )}. (3.10) 

Since a 2 and (03+05) are minimum at Nfr ~ 3 (see Table I), the decay is obviously dominated 
by the electroweak penguin transition when N° S (LL) = N° S (LR) = 3. Numerically, we find 
at the amplitude level 

tree : QCD penguin : electroweak penguin = 0.28 : 1 : —2.72 . (3-H) 
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It is clear that although Rw = 0.79 for N° s = 3, the decays B s — > ujt( and B s — > uxfi are 
actually dominated by the electroweak penguin. 

The branching ratios for the class-V and -VI modes shown in (3.9) depend strongly on the 
value of N% . As pointed out in Sec. II, the preferred values for the effective number of colors 
are N^ S (LL) « 2 and N° S (LR) ~ 5. We believe that the former will be confirmed soon by 
the forthcoming measurements of B — > nn, np, ■ ■ ■. However, the branching ratios for some 
of the decay modes, e.g. B s — > ur], urj', (f)rj, become very small at the values of Nf* given by 



Eq. (2.13). As suggested in [22], these decays involve large cancellation among competing 
amplitudes and they may receive significant contributions from annihilation and/or final- 
state interactions. 

As noted in passing, class-IV modes involve the QCD penguin parameters a 4 and a 6 in 
the combination a 4 + Ra 6 , where R > for B s — > P a Pb, R = for P a Vb and V a Vb final states, 
and R < for B s — > V a Pb, where Pb or V& is factorizable under the factorization assumption. 
Therefore, the decay rates of class-IV decays are expected to follow the pattern: 

T(B S - P a P b ) > T(B S - P a V b ) ~ T(B S - V a V b ) > T(B S ^ V a P b ), (3.12) 

as a consequence of various possibilities of interference between the penguin terms charac- 
terized by the effective coefficients a 4 and a 6 . From Tables III-V, we see that 

T(B S -> K + K~) > T(B S -> K + K*-) > T(B S -> K*+K*~) > T(B S -> K + *R-), 

t(b s -> ir°Z°) > r(5 s -> k°k*°) > r(s s -> ir*°K*°) > r(s s -> ir*°K°). (3.13) 

Note that the pattern r(S — > P a V&) > T(i? — > P a Pb), which is often seen in tree-dominated 
decays, for example r(5 s — > K + p~) > T(B S — > K + tt~), occurs because of the larger spin 
phase space available to the former due to the existence of three different polarization states 
for the vector meson. On the contrary, the hierarchy (3.13) implies that the spin phase-space 
suppression of the penguin-dominated decay B s — > P a Pb over B s — ► P a Vb or B s — > V a P b 
is overcome by the constructive interference between penguin amplitudes in the former. 
Recall that the coefficient R is obtained by applying equations of motion to the hadronic 
matrix elements of pseudoscalar densities induced by penguin operators. Hence, a test of 
the hierarchy shown in (3.13) is important for understanding the calculation of the penguin 
matrix element. | 

Among the 39 charmless two-body decay modes of the B s meson, we find that only seven 
of them have branching ratios at the level of 1CT 5 : 

B S ^K + K~, K°K°, TjTj 1 , r/r/, K + p~, K + *p~, <jxf>. (3.14) 

It is interesting to note that among the two-body rare decays of B~ and Bj, the class- VI 



decays B — > r{K and B<i — > rj'K have the largest branching ratios [pqj : 



§For a direct estimate of R using the perturbative QCD method rather than the equation of 
motion, see 
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B{B± -> r/K*) = [6.5t\i ± 0.9J x 10~ 5 , 

B{B d -> t/K ) = (4.7^.7 ± g) x 10 -5 ( 315 ) 



The decay rate of B~ — > r(K~ and P^ — > r]'K° is large because they receive two different sets 
of penguin contributions proportional to a 4 + Ra 6 with R > 0. By contrast, VP, VV modes 
in charm decays or bottom decays involving charmed mesons usually have larger branching 
ratios than the PP mode. Because of the strange quark content of the B s , one will expect 
that the decay B s — > i]r( or B s — > 7/7/, the B s counterpart of B d — > rj'K , is the dominant 
two-body P s decay. Our calculation indicates that while the branching ratio of P s — > 7777' is 
large, 

/3(P S -> 777/) «2x 1(T 5 for Nf(LL) = 2, Nf(LR) = 5, (3.16) 

it is only slightly larger than that of other decay modes listed in (3.14), see Tables III-V. 

What is the role played by the intrinsic charm content of the rf to the hadronic charmless 
B s decay ? Just as the case of B — ► rj'K, B s — > 77^7/' receives an internal TV-emission 
contribution coming from the Cabibbo-allowed process b — > ccs followed by a conversion 
of the cc pair into the 7/ via gluon exchanges. Although the charm content of the 77' is 
a priori expected to be small, its contribution is potentially important because the CKM 
mixing angle Vf>V c * is of the same order of magnitude as that of the penguin amplitude [cf. 
Eqs. (A10,A11)] and yet its effective coefficient a 2 is larger than the penguin coefficients by 
an order of magnitude. Since a 2 depends strongly on N° S (LL) (see Table I), the contribution 
of cc — > t]' is sensitive to the variation of N° (LL). It is easy to check that the 77' charm 
content contributes in the same direction as the penguin terms at 1/N^ S (LL) > 0.28 where 
a 2 > 0, while it contributes destructively at 1 / 7V ( ? ff (LL) < 0.28 where a 2 becomes negative. 
In order to explain the abnormally large branching ratio of B —>■ r]'K, an enhancement from 
the cc — > t]' mechanism is certainly welcome in order to improve the discrepancy between 
theory and experiment. This provides another strong support for N^ S (LL) w 2. Note that 
a similar mechanism explains the recent measurement of B~ — > r] c K~ ||37|| . 

It turns out that the effect of the cc admixture in the 77' is more important for B s — > rj'rj' 
than for B s — > tjt]'. It is clear from Eq. (Al) that the destructive interference between 
x {B s r hV >) x fcpBsv and x {B sV ',v) K f^F^' in the decay amplitude of B s -> 7777', recalling 

that the form factors Fq sV and F sV have opposite signs, renders the contribution of cc — > v[ 
smaller for B s — > r)T) ; . 

A very recent CLEO reanalysis of B — > rj'K using a data sample 80% larger than in 



previous studies yields the preliminary results p8 



B{B ± -> rfK*) = (jAtH ± l.O) x 10" 5 , 

B(B d -> 77'ir ) = (5.9±i;| ± 0.9) x 10" 5 , (3.17) 

suggesting that the original measurements (3.15) were not an upward statistical fluctuation. 
This result certainly favors a slightly larger f c (l) in magnitude than that used in (2.27). 



In fact, a more sophisticated theoretical calculation gives fZ = —(12.3 ~ 18.4) MeV [p9| , 
which is consistent with all the known phenomenological constraints. This value of f£, will 
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lead to an enhanced decay rate for B — > rj'K. Numerically, we find that for N^ S (LL) = 2, 
Nf(LR) = 5 and = -15 MeV, 

J3(B S -> 7777') = 2.2 x 10~ 5 , £(5 a -»■ 7/7/) = 1-8 x 10~ 5 , (3.18) 

to be compared with 

/3(v3 s -> 7777') = 1.8 x 1(T 5 , B{B S -> 7/7/) = 1.2 x 10~ 5 , (3.19) 

in the absence of the intrinsic charm content of the 7/. 

Finally, we should point out the uncertainties associated with our predictions. Thus far, 
we have neglected H^-annihilation, space-like penguin diagrams, and final-state interactions; 
all of them are difficult to estimate. It is argued in P2| that these effects may play an essential 
role for our class-V and -VI decay modes. Other major sources of uncertainties come from the 
form factors and their q 2 dependence, the running quark masses at the scale m&, the virtual 
gluon's momentum in the penguin diagram, and the values for the Wolfenstein parameters 
p and 77. 

IV. CONCLUSIONS 

Using the next-to-leading oder QCD-corrected effective Hamiltonian, we have systemat- 
ically studied hadronic charmless two-body decays of B s mesons within the framework of 
generalized factorization. Nonfactorizable effects are parametrized in terms of N^ G (LL) and 
N° G (LR), the effective numbers of colors arising from (V — A)(V — A) and (V — A)(V + A) 4- 
quark operators, respectively. The branching ratios are calculated as a function of N^ S (LR) 
with two different considerations for N^ S (LL): (i) N° S (LL) being fixed at the value of 2, 
and (ii) N^ S (LL) = N^ S (LR). Depending on the sensitivity of the effective coefficients af 
on N% , we have classified the tree and penguin transitions into six different classes. Our 
results are: 

1. The decays B s — > r]7c, i]'tt, rjp, rj'p, (pir, <pp receive contributions only from the tree 
and electroweak penguin diagrams and are completely dominated by the latter. A 
measurement of them can be utilized to fix the effective electroweak penguin parameter 
ag. For N^ S (LL) = 2, we found that electroweak penguin contributions account for 
85% of their decay rates. Their branching ratios, though small [in the range of (0.4 — 
4.0) x 10~ 7 ], could be accessible at hadron colliders with large b production. 

2. For class-V and -VI penguin-dominated modes: B s — > ur/, urf, (pi], ur/, K<f), K*<f), (ficj), 
electroweak penguin corrections, depending strongly on N* , are as significant as QCD 
penguin effects and can even play a dominant role. 

3. Current experimental information on B~ — > ujtt~ and B° — > 7r + 7r - favors a small 
N% (LL), that is, N^ S (LL) w 2, which is also consistent with the nonfactorizable term 
extracted from B — > (D, D*)(tt, p) decays, N^ S (B — > Dir) rs 2. We have argued that 
the preferred value for the effective number of colors N° S (LR) is N^ S (LR) ~ 5. 
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4. Because of various possibilities of interference between the penguin amplitudes gov- 
erned by the QCD penguin parameters and a^, the decay rates of class-IV decays 
follow the pattern: T(B S -> P a P b ) > T(B S -> P a V b ) ~ T(5 S -> KH) > r(P s -> V P 6 ), 
where P a = K + , P b = K~ or P a = K°, P b = K . A test of this hierarchy is important 
to probe the penguin mechanism. 

5. The decay B — > rj'K is known to have the largest branching ratios in the two-body 
hadronic charmless B~~ and Bd decays. Its analogue in the B s system, namely B s — > r/r/' 
has a branching ratio of order 2 x 10~ 5 , but it is only slightly larger than that of 
rj'rj', K* + p~ , K + K~ , K°K decay modes, which have the branching ratios of order 10~ 5 . 

6. The recent CLEO reanalysis of B — > rj'K favors a slightly large decay constant /£. Us- 
ing f^i = — 15 MeV, which is consistent with all the known theoretical and phenomeno- 
logical constraints, we found that the intrinsic charm content of the rf is important for 
B s — > rj'rj 1 ', but less significant for B s — > 7777'. 
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APPENDIX 
A. The B s — > PP decay amplitudes 

For B s — > PP decays, we use X^ BsPl,p ^ to denote the factorizable amplitude with the 
meson P 2 being factored out. Explicitly, 

x (b 3 p u p 2) = ( p 2 |(^3 )v _jo)(P 1 |( ? - 1 6) v _jP s ) = i/ ft (ro|. (Al) 

For a neutral Pi with the quark content N(qq + •••), where N is a normalization constant, 

x (M,flO = (P 2 |(^)^_jO)(P 1 |( g - 1 6) v _jP s ) = i/4(m| s - m^J^K). (A2) 

As an example, the factorizable amplitudes X^ BsV ' ,K ^ and X^ BaK,v '^ of the decay P s — > if ?/ 
read 

X (iW,JO = (Z°|( S -d) v _jO)( ?7 '|(J6) v _jP s ) = ?/*(< -m^o^'Kc), 

= ( 77 '|(^),_JO)(F°|( S -6) v _ A |P s ) = - ( A3 ) 

For simplicity, VF-annihilation, space-like penguins and final-state interactions are not in- 
cluded in the decay amplitudes given below. 

1. b — > d processes: 
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A(B S -> K+n-) = ^\v uh V: d a x - V tb V t * d 



a 4 + a w + 



+ 2(a 6 + a 8 ) 



(m„ + m d )(m b - m u ) 



]^ X (B s K+,n-)^ 



(A4) 



3. , 1 

a 4 + 2^ _fl7 + fl9 ) + 2 ai0 



- 2(a 6 - -a 8 ) 



2 (m d + m d )(m b -m d )_ 



>X^°\ (A5) 



- W 



K 



1 1 

a 4 -aio + 2(a 6 -a 8 ) ^ + ^ ^ _ ^ . 



+ (a s - a 5 - a 7 + a 9 ) xf^ ) + (a 3 - a 5 + \a 7 - ^a 9 )Xf 



+ (a 3 -a 5 - a 7 + a 9 )Xf- 



( 111 

+ a 3 + a 4 -a 5 + -a 7 - -a 9 - -ai 



1 



+ («6 - ^ a s) 



2 m s (m b — m s 



x {b,k, v <!)) 



(A6) 



2. 6 — > s processes: 



A(B S 



K + K~) = ^ ] V ub V: s a 1 -V th V; s 



a 4 + aio + 



+ 2(a 6 + a 8 )- 



A' 



m u + m s )(m b - m u ) 



>X (B S K+,K-)^ 



(AT) 



A(B S -> A°) = ^§{*W, 02 - W s 



;(-a 7 + a 9 ) 



| X (B S ,('),.0) ; (Ag) 



where 



(A9) 
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- v tb v t 



ts 



' 111 

a 3 + a 4 - a 5 + -a 7 - -a 9 - -ai 



1 



+ (a 6 - -a & ) 



rat 



2 m s (m 6 - m s ) \ /* 



+ (2a 3 - 2a 5 - ^a 7 + ^a 9 ) xf*™') + (a 3 - o 5 - a 7 + a 9 )Xj B '™'> 



/ 111 

+ a 3 + a 4 -a 5 + -a 7 - -a 9 - -ai 



+ K - ^as) 



mz 



2~ oy m s (m 6 -m s ) ^ /, s jj s 



+ ^2a 3 - 2a 5 - ^a 7 + ^a 9 ) jtfW") + (a 3 - a 5 - a 7 + a 9 )X c ^'^ |, 



(A10) 



A(P S - rfrf) = ^2lv ub V: s a 2 X^'^ + V cb V c * s a 2 X^'^ 



V tb V t 



ts 



1 1 1 

a 3 + a 4 - a 5 + -a 7 - -a g - -ai 



+ (a 6 - -a 8 )- 



2 m s (m b + m s ) \ 
+ (2a 3 - 2a 5 - l -a 7 + ^a 9 ) ^-"') + (a 3 - a 5 - a 7 + a 9 )Xf ^-"') j. (All) 

The amplitude of B s — > 7777 is obtained from A(P S — > 7/77') by the replacement 7/ — > 77. 
3. pure penguin process: 



A{B S 



K°K°) = ^i-v tb v; s 



1 

a 4 — 2 0,10 



> 2 K 



1 777." 



>X (B S K°,K°)_ 



(A12) 



B. The P s — > VP decay amplitudes 

The factorizable amplitudes of P s — > VP decays have the form: 

x (b s p, V) = (V|(g 2 g 3 ) v _jO)(P|(g 1 6) v _j5 s ) = 2/ v m y F^ P (m 2 y )(5 • Pfl J, 
x (b s v,p) = (P|(g 2g3 ) v _J )(\/|(g 1 6) v ,_jP s ) = 2/ pm ^rK)( £ .j,J. 



(Bl) 



For example, the factorizable terms X^ BaV '' K *^ and X^ BsK *' v,s> of P s — > K*rj' decay are given 
by 
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X (W*) = (K^\(3u) v _ A \0)(rf\(vb) v _ A \B a ) = 2f K .m K .F 1 B "<(m] c .)(e-p B .), 
X (B S K*, V ') = {r f mv _j Q){K ^ mv _ A \B s ) = 2f«m K *A*° K \ml,){s-p Bs ). (B2) 



I. b — > d processes: 



2(a 6 + a 8 ) 



2 



Cl 4 + Q10 



ml 



[m u + m d )(m b + m u ) 



>x (b s k+*,k-) 



(B3) 



A(B S -> = ^jWd ai - ^(04 + am) }* (fl ' 



K+,p~) 



(B4) 



+ 2(a 6 2 a8 \m d {m h + m d ) 



3 3 1 

>x (B s i^>V°) 



(B5) 



- K°p°) = ^{v ub V: d a 2 - V tb V; d ( - a 4 + |a 7 + ^a 9 + i fll0 ) }xf .*V>, 



(B6) 



- K°lo) = ^V ub V: d a 2 - V tb V* d (2a 3 + a 4 + 2a 5 



(B7) 



A{B S 



- W 



(a 4 - iaio)X (Bsr? + fa 3 + a 4 - a 5 + ^a 7 - ^a 9 - -a 10 



- K - ~a 8 ) 



2 
v 



'JL - 1 I r , I X (Bs ^* ,,? ) 
2^m s (m fe + m s ) IP, I " ' 



+ (a 3 - a 5 - a 7 + a 9 ) 
+ (a 3 -a 5 -aT + a 9 )lf r/>) |. 



^ + (a 3 - a 5 + ia 7 - k)^'^) 



(B8) 
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2. b — > s processes: 



A(5 S - x+k-*) = ^| v ub v: s 0l - v^;(a 4 + a 10 ) J*^.*-), 



V2 
2(a 6 + a 8 ) 



a 4 + aio 



m 



A' 



(m s + m u )(m b + m u ) 



]^ X (B S K+*,K-)^ 



a{b s -> A ) = ^|k 6 C «2 - w 



-(a 7 + a 9 ) 



2(a 3 + a 5 ) + -(a 7 + a 9 ) 



>Xi B ° J > ( ' ) >"\ 



where 



-(-07 + 09) 



>xi B >*>*°\ 



A(B S - H n ) = ^[vu b V: s a 2 Xi B ^ + V cb V: s a 2 X^^ 



(a 6 



1 1 1 

a 3 + «4 - 05 + ^ a 7 _ ~ 2 fll ° 



1 x Ko 



fu 

1_ Vl] ) X (Bs4>V>) 



2 m s (m b + m s ) y f' {l) 
+ (2a 3 - 2a 5 - ia 7 + ^a 9 ) X< B '*>" (,) > + (a 3 - a 5 - a 7 + a 9 )Xf » 

+ (a 3 + a 4 + a 5 - ^a 7 - ^9 - ^a^X^^ |. 



3. pure penguin processes: 



-^ Ws (a 4 -I ai0 )x^ 



K°,K *) 
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A(B S -> K°*K 



G_f_ 
V2 



-^=v tb v* 



1 



a 4 — -aio 



2 («6 2 a§ ' ) (m s + m d )(m 6 + m d ) 



V2 



a 3 + a 5 - -(a 7 + a 9 ) 



+ 



m 



A' 



1 1 X 

«4 - -«io - 2(a 6 - -a 8 )- ■ 

2 2 (m s + m d )(m b + m d )_ 



C. The B s — > VV decay amplitudes 

The factorizable amplitude of 5 S — > W decays has the form: 



X (B s VuV.) = lfv2mv2 



e\ ■ e* 2 )(m Bs + m Vl )A^(m 2 V2 ) 



i _ ^ rt _ ^ // fv t) 



(m Bs + m Vl ) 



{m Bs +m Vl ) 



I. b 



d processes: 
A(B. -> K + *p-) 



ax - ^(04 + a 10 ) jx^^"), 



CrF -iv ub v: d a 2 -v tb v: d { 



V2\ 



3 3 
a 4 + -07 + + o 



A(£? s - K°*u) = ^V uh v: d a 2 - V tb V t * d (2a 3 + a 4 + 2a 5 

+ ^ 7 +^9-^xo)}x^-). 



2. b — > s processes: 



Op 
V2' 



V tb V*(a, + a w )}x^ K+ *> K -*\ 
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A(B S -> p°0) = ^| V ub V: s a 2 - V tb V* 



(a 7 + a 9 ) 



>x (B s <j>,pO) 



U ) 



A(B S -> w = ^|K 6 C «2 - V»V t 



2(a 3 + a 5 ) + -(a 7 + a 9 ) 



> J£(B B <t>,UJ 



pure penguin processes: 



A(5 S - K°*Z *) 



_ 7f ^ 2 ^ 



+ (a 4 -±a 10 )x< fl -^>}, 



^3 + ^5 - + °g) 



a 3 + a 4 + a 5 - -(a 7 + a 9 + ai ) 
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Table III. Branching ratios (in units of 10 6 ) averaged over CP-conjugate modes for charmless 
B s -»• PP decays. Predictions are for k 2 = ml/2, 77 = 0.34, p = 0.16, and Nf(LR) = 2,3,5,oo 
with N^ S (LL) being fixed to be 2 in the first case and treated to be the same as N^ S (LR) in the 
second case. We use the BSW model for form factors [see (2.28)]. 



Decay 


Class 




Nf(LL) = 2 




Nf(LL) 


= Nf(LR) 


2 


3 


5 


00 


2 


3 


5 


00 


B s -► K+7T- 


I 


6.64 


6.66 


6.67 


6.70 


6.64 


7.38 


8.01 


8.99 


B s - K°ir° 


II 


0.24 


0.24 


0.25 


0.25 


0.24 


0.08 


0.12 


0.46 


B s -► K + K~ 


IV 


9.88 


10.9 


10.9 


11.6 


9.88 


10.9 


11.7 


12.9 


B s - 


IV 


10.3 


10.9 


11.4 


12.1 


10.3 


12.0 


13.5 


15.8 


5, A' 


V 


0.04 


0.04 


0.04 


0.04 


0.04 


0.05 


0.06 


0.09 


5 S 7r°r/ 


V 


0.04 


0.04 


0.04 


0.04 


0.04 


0.05 


0.06 


0.09 




VI 


0.63 


0.86 


1.06 


1.42 


0.63 


0.54 


0.57 


0.76 


U s -► K°ry 


VI 


0.81 


0.84 


0.87 


0.91 


0.81 


0.82 


0.96 


1.39 


-B s -> W 


VI 


12.5 


16.3 


19.6 


25.3 


12.5 


14.4 


15.9 


18.5 


5 S 7/7/ 


VI 


6.28 


10.3 


14.3 


21.4 


6.28 


6.80 


7.23 


7.91 


5 S -> 777/ 


VI 


5.30 


4.80 


4.41 


3.89 


5.30 


6.23 


7.05 


8.37 
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Table IV. Branching ratios (in units of 10 6 ) averaged over CP-conjugate modes for charmless 
B s -» VP decays. Predictions are for k 2 = ml/2, rj = 0.34, p = 0.16, and Nf(LR) = 2,3,5,oo 
with N° S (LL) being fixed to be 2 in the first case and treated to be the same as N^ S (LR) in 
the second case. For decay modes involving the B s — * K* or B s — > <j) transition, we use two 
different models for form factors: the BSW model [jl8| (the upper entry) and the light-cone sum 
rule approach [34| (the lower entry). 



Decay 


Class 




Nf(LL) 


= 2 






Nf(LL) = 


Nf(LR) 




2 


3 


5 


oo 


2 


3 


5 


oo 


B s ^ 


K* + ir- 


I 


4.30 


4.30 


4.30 


4.30 


4.30 


4.79 


5.20 


5.84 








4.98 


4.98 


4.98 


4.98 


4.98 


5.55 


6.02 


6.76 


B s ^ 


K+p- 


I 


17.2 


17.2 


17.2 


17.2 


17.2 


19.2 


20.8 


23.4 


B s ^ 


K°*7T 


II 


0.14 


0.14 


0.14 


0.15 


0.14 


0.01 


0.02 


0.23 








0.17 


0.17 


0.17 


0.17 


0.17 


0.01 


0.02 


0.27 


B s ^ 


K°p° 


II 


0.55 


0.55 


0.55 


0.55 


0.55 


0.06 


0.13 


1.00 


B s ^ 




II, VI 


0.10 


0.14 


0.18 


0.26 


0.10 


0.04 


0.04 


0.15 








0.11 


0.15 


0.20 


0.29 


0.11 


0.03 


0.04 


0.17 


B s ^ 


K*° V 


II,VI 


0.18 


0.18 


0.18 


0.17 


0.18 


0.13 


0.17 


0.39 








0.20 


0.19 


0.19 


0.19 


0.20 


0.13 


0.18 


0.42 


B s ^ 


K°u 


II, VI 


0.71 


0.60 


0.53 


0.46 


0.71 


0.11 


0.07 


0.77 


B s ^ 


K+*K~ 


IV 


0.68 


0.78 


0.87 


1.01 


0.68 


0.75 


0.80 


0.88 








0.79 


0.90 


1.00 


1.16 


0.79 


0.86 


0.92 


1.02 


B s ^ 


K°*K° 


IV 


0.26 


0.34 


0.41 


0.53 


0.26 


0.20 


0.15 


0.10 








0.31 


0.40 


0.48 


0.62 


0.31 


0.23 


0.18 


0.11 


B s ^ 


K+K-* 


IV 


3.40 


3.40 


3.40 


3.56 


3.40 


3.77 


4.07 


4.55 


B s -> 


K°K°* 


IV 


3.28 


3.28 


3.28 


3.28 


3.28 


4.15 


4.92 


6.21 


B s ^ 


vrV 


V 


0.18 


0.18 


0.18 


0.17 


0.18 


0.22 


0.27 


0.40 








0.35 


0.34 


0.34 


0.33 


0.35 


0.42 


0.53 


0.78 


B s ^ 


prj' 


V 


0.11 


0.11 


0.11 


0.11 


0.11 


0.13 


0.17 


0.26 


B s ^ 


m 


V 


0.11 


0.11 


0.11 


0.11 


0.11 


0.14 


0.18 


0.26 


B s ^ 


LOT]' 


V 


0.79 


0.18 


0.01 


0.31 


0.79 


0.004 


0.36 


2.52 


B s ^ 




V 


0.80 


0.18 


0.01 


0.31 


0.80 


0.004 


0.36 


2.56 


B s ^ 


H 


VI 


1.06 


1.18 


1.28 


1.45 


1.06 


0.27 


0.22 


1.11 








0.55 


0.86 


1.20 


1.86 


0.55 


0.31 


0.75 


2.45 


B s ^ 


<t>V 


VI 


2.03 


0.79 


0.25 


0.20 


2.03 


0.91 


0.34 


0.04 








1.43 


0.41 


0.15 


0.69 


1.43 


0.58 


0.19 


0.09 


B s ^ 


K°cp 


VI 


0.002 


0.01 


0.04 


0.13 


0.002 


0.03 


0.10 


0.29 








0.004 


0.03 


0.07 


0.19 


0.004 


0.04 


0.12 


0.32 
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Table V. Same as Table IV except for B s — > VV decays. 



Decay 


Class 




Nf(LL) 


= 2 






Nf(LL) = 


Nf(LR) 




2 


3 


5 


oo 


2 


3 


5 


oo 


B s K+*p~ 


I 


12.5 


12.5 


12.5 


12.5 


12.5 


13.9 


15.0 


16.9 






14.4 


14.4 


14.4 


14.4 


14.4 


16.0 


17.4 


19.5 


B s -► K°*p° 


II 


0.40 


0.40 


0.40 


0.40 


0.40 


0.044 


0.094 


0.72 






0.46 


0.46 


0.46 


0.46 


0.46 


0.051 


0.11 


0.84 


B s -► 


II,VI 


0.26 


0.21 


0.19 


0.17 


0.26 


0.04 


0.02 


0.28 






0.30 


0.25 


0.22 


0.19 


0.30 


0.044 


0.031 


0.32 


5 S -► K+*K-* 


IV 


2.53 


2.53 


2.53 


2.53 


2.53 


2.80 


3.03 


3.38 






2.91 


2.91 


2.91 


2.91 


2.91 


3.22 


3.48 


3.88 


B s -► 


IV 


2.44 


2.44 


2.44 


2.44 


2.44 


3.09 


3.66 


4.62 






2.80 


2.80 


2.80 


2.80 


2.80 


3.55 


4.21 


5.30 


s s - 


V 


0.17 


0.18 


0.18 


0.18 


0.17 


0.22 


0.28 


0.41 






0.33 


0.34 


0.34 


0.35 


0.33 


0.42 


0.53 


0.79 


-B s — > 


V 


0.65 


0.15 


0.01 


0.25 


0.65 


0.004 


0.30 


2.09 






1.22 


0.27 


0.02 


0.48 


1.22 


0.007 


0.56 


3.92 


5 S -► 


VI 


0.007 


0.049 


0.10 


0.19 


0.007 


0.13 


0.28 


0.57 






0.014 


0.098 


0.17 


0.30 


0.014 


0.22 


0.43 


0.86 


S s -» 00 


VI 


13.8 


8.77 


5.57 


2.15 


13.8 


7.15 


3.40 


0.37 






25.1 


15.9 


10.1 


3.91 


25.1 


13.0 


6.18 


0.68 
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Table VI. Fractions of non-electroweak penguin contributions to the branching ratios of penguin- 
dominated two-body B s decays, as defined by Eq. (3.8). Predictions are for k 2 = ml/2, rj = 
0.34, p = 0.16, and Nf(LR) = 2,3,5,oo with Nf(LL) being fixed to be 2 in the first case and 
treated to be the same as N° S (LR) in the second case. We use the BSW model for form factors. 



Decay 




l\ c [LL) 


— Z 






iV c (LL) 


— iV c (LK) 




z 


Q 

o 


E 




oo 


Z 


q 


c 



OO 


~d / 


0.15 


0.15 


0.16 


0.16 


0.15 


0.007 


0.01 


/\ -1-1 

0.11 


B s 7T°7/ 


0.15 


0.15 


0.16 


0.16 


0.15 


0.007 


0.01 


0.11 


B s ttV 


0.15 


0.15 


0.16 


0.16 


0.15 


0.007 


0.01 


0.11 


B s A' 


0.15 


0.15 


0.15 


0.14 


0.15 


0.007 


0.01 


0.11 


5 S A 


0.15 


0.15 


0.15 


0.14 


0.15 


0.007 


0.01 


0.11 


5 S -► wr/ 


0.78 


0.57 


1.63 


1.43 


0.78 


0.79 


1.42 


1.16 


£ s -> cjt/ 


0.78 


0.57 


1.63 


1.43 


0.78 


0.79 


1.42 


1.16 


5 S -► 0r/ 


1.73 


1.70 


1.69 


1.65 


1.73 


1.93 


0.63 


0.61 


5 S -> <M 


1.71 


2.16 


3.01 


0.39 


1.71 


2.00 


2.58 


2.81 


5 S - K> 


3.25 


0.23 


0.49 


0.07 


3.25 


0.43 


0.68 


0.82 


5 S -»• K* + K*- 


0.87 


0.87 


0.87 


0.87 


0.87 


0.94 


1.00 


1.08 


B s - /r°iT° 


1.08 


1.08 


1.08 


1.08 


1.08 


1.03 


1.00 


0.96 


s s - A 


0.15 


0.14 


0.14 


0.14 


0.15 


0.006 


0.01 


0.11 


-B s — > u;0 


0.78 


0.57 


1.63 


1.43 


0.78 


0.79 


1.42 


1.16 


5 S -► 


3.82 


0.55 


0.76 


0.86 


3.82 


0.75 


0.84 


0.87 


5 S -► 00 


1.25 


1.32 


1.41 


1.69 


1.25 


1.32 


1.43 


2.19 
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